Droplet traffic in microfluidic networks: 
A simple model for understanding and designing 
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We propose a simple model to analyze the traffic of droplets in microfluidic "dual networks" . Such 
functional networks which consist of two types of channels, namely those accessible or forbidden to 
droplets, often display a complex behavior characteristic of dynamical systems. By focusing on three 
recently proposed configurations, we offer an explanation for their remarkable behavior. Addition- 
ally, the model allows us to predict the behavior in different parameter regimes. A verification will 
clarify fundamental issues, such as the network symmetry, the role of the driving conditions, and of 
the occurrence of reversible behavior. The model lends itself to a fast numerical implementation, 
thus can help designing devices, identifying parameter windows where the behavior is sufficiently 
robust for a devices to be practically useful, and exploring new functionalities. 



Injecting immiscible fluids into a microfluidic sys- 
tem has been found to permit the controlled generation 
of droplets or bubbles [l|. Rapidly, droplet micro flu- 
idics has developed beyond early observations because 
of its promises for a diversity of fields @, Q: studies 
of chemical kinetics and time-controlled synthesis [3], 
high-throughput screening fabrication of remarkable 
colloidal objects combinatorial chemistry and 
even non-electronic coding or computing functions Q • 
These tasks require mastering the generation of droplets, 
their traffic in potentially complex networks, and some- 
times their splitting or merging. One can either employ 
passive means by relying on hydrodynamics and capillar- 
ity alone 0, 0, [j, 0, 0] , or integrate active actuators (elec- 
tric, dielectric, pneumatic) in the devices [2,l9(. Difficul- 
ties often arise even in simple geometries if the complex 
behavior of dynamical systems (period doubling, chaotic 
behavior, etc.) occurs |Ql, 1 1 0f| - 

We here focus on the traffic of droplets only, with- 
out splitting or merging, in what we name dual net- 
works, i. e. networks combining channels that can be ac- 
cessed by droplets (transport channels) and others that 
cannot ("bypass" channels), see Fig. [TJ This family 
of networks encompasses many recently proposed solu- 
tions @, EE El- Building on recent works [111, El, 
we introduce a simple numerical tool to analyze the col- 
lective motion of droplets in such a dual network. We 
show this tool to be powerful in both illuminating fun- 
damentals ruling the traffic (including considerations of 
symmetry and reversibility) and as a means to predict the 
capability and robustness of a given design in performing 
a given function. In particular, we illustrate these ideas 
by analyzing three recently proposed devices as displayed 
in Fig. [TJ a reversible encoding loop [8( (jTJa), a ladder for 
synchronizing droplets between parallel channels ([It 1 ), 
and a bypassed T-junction for symmetric splitting of a 
droplet train (H |JTfc). 

Simple Model - We model both the flow of the contin- 
uous phase and the traffic of droplets in an arbitrary dual 
network of transport channels and bypasses, represented 
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FIG. 1: Graph representations of three recently proposed dual 
networks: (a) a loop Q; (b) a ladder @|; (c) a bypassed T- 
junction [12]. Full and dotted lines correspond to transport 
channels (accessible to droplets) and "bypass" channels (for- 
bidden to droplets). Droplets are indicated by circles. 



as a simple directed graph of nodes and links as in Fig.[TJ 
The instantaneous flow state is given by a pressure P a at 
each node and a total flowrate Qi in each channel. These 
variables are related by equations similar to Kirchhoff's 
rules for electrical circuits. Mass conservation at node 
a, taking into account a possible external injection Qgxt) 
reads 
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where e; = ±1 account for the orientation of the link. 
In channel i the flowrate Qi is related to the pressure 
difference between its two end nodes QWj) and a ut(j). 
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Equation ([2]) defines the hydrodynamic resistance Ri of 
channel i. For flows without droplets at low Reynolds 
number it reduces to a constant value Ri, proportional 
to the viscosity and to the channel length Li. Ri oc Li. 
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Droplets or bubbles in the channel significantly modify 
the picture in a quite complex way [3], OB]!. A simplifi- 
cation occurs if the droplets are sufficiently distant from 
one another (typically a few channel diameters away) so 
that their flow perturbations do not interact. In the men- 
tioned experiments the ratio of distance to diameter is ap- 
prox. 10 (loop and ladder) and 5 (bypass junction). Each 
droplet then yields an additive increment Rd to the resis- 
tance. In this paper, we make the further simplification 
that all droplets are equal and so are the cross-sections S 
of the channels transporting droplets (these restrictions 
can be easily lifted in the numerical model), so that for 
Hi droplets in channel i, 



Ri — Ri 
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Identical droplets with no interaction move all at the 
same velocity Vi in a channel with a given flow rate Qi, 
which we write with a proportionality factor (3/ S, 

Vi = 0Qi/S. (4) 

In general, both Rd and /3 will vary with Qi [l3|j [T3 |. 
The simplifications in ^ and Q, however, have proved 
to hold in experimental realizations [llj |. 

The flow state, consisting of all pressure and flowrate 
values, is uniquely determined by the set of linear equa- 
tions JTJ) and @ when complemented by appropriate 
boundary conditions at the end-nodes, i. e. specification 
of either the pressure or the externally injected flowrate. 
An important point is that this flow state persists un- 
til a droplet arrives at a node and enters another chan- 
nel, thereby modifying the set of resistances. To de- 
scribe droplet traffic, we thus track the locus of the 
droplets (modeled as points) along the channels. When 
a droplet arrives at a node where there are several chan- 
nels to choose from, a selection rule is required. We fo- 
cus here on locally symmetric T-junctions for which the 
droplets enter the channel with the instantaneous larger 
flow rate [UGl. 

Altogether, we propose the following algorithm: 
(i) Compute the pressures and flowrates for a given set 
of droplet positions using Eqs. (JTJ) and (J5J); (ii) determine 
the droplet velocities according to Eq. (in) determine 
the next time when either a droplet arrives at a junction 
or a new one is injected, and advance all droplets to this 
time; (iv) decide which route the triggering droplet takes, 
update the channel resistances of the affected channels, 
and return to step (i). 

Clearly, step (iii) is the big time saver that allows fast 
computations of rather involved histories. This scheme 
can be straightforwardly adapted to networks with chan- 
nels of different cross-sections and to other junction ge- 
ometries, but the present version is sufficient to revisit re- 
cently proposed devices and discuss fundamental issues. 

Reversibility - Fuerstman et al. presented a device 
that allows reversible "encoding" of information in in- 
terdroplet distances: Two long channels, used to store 



droplet series, are connected by a slightly asymmetric 
loop (Fig. QJi). A regular train of droplets, injected 
through one channel into the loop, exits the loop as a 
complex modulated sequence. The interdroplet distances 
result from the "decisions" taken by the droplets when 
entering either branch of the loop. Each decision is af- 
fected by the presence of previous droplets in the loop. 
The slight asymmetry of the loop gives rise to periods of 
2, 3, ... 7 in the outgoing sequence, depending on the ini- 
tial interdroplet distance Q . When we apply our simple 
model using parameters reasonably related to those of 
the experimental realization, we do obtain a comparable 
cascade of periods (2,3,5,7,9 etc.), as a characteristics 
of dynamical systems at work in this network [l5l ]. 

A remarkable result by Fuerstman etal. [8| is that 
upon reversal of the flow direction, the starting sequence 
of droplets is restored, i. e. the system exhibits reversibil- 
ity. This is far from being obvious, keeping in mind that 
the selection rule is not intrinsically reversible. 

To investigate the question of reversibility further, in 
our numerics we send a periodic train of droplets through 
the loop until N droplets have reached the other side, 
then invert the flow direction and send them back. The 
resulting (N — 1) interdroplet distances after this return 
trip through the loop are plotted in Fig. for sev- 
eral values of the initial distance A, ranging from to 
the maximum for which two droplets interact, A max = 
Le[l + (Re + Rd)/R u ] (see Fig. QJ, for the parameters). 
The most striking feature of Fig. [2^ is that it predicts 
windows of values of A for which the system behaves re- 
versibly and other windows for which it does not (points 
not on the diagonal). It also sheds light on the impor- 




FIG. 2: Final distances Ax n between 70 droplets after a re- 
turn journey through the loop of the device in Fig. [jj,, as 
a function of the initial period A. Reversibility is observed 
where all 69 symbols collapse on the diagonal, see also the 
zoomed panel b. The results differ for a finite droplet train 
(panels c and d), where also the number of droplets matters. 
The parameters L u /Le = 1.112, Rd/ Re = 1.5, are compatible 
with data from Ref. [3 [Hj]. 
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tance of the precise conditions at flow inversion: For 
the same parameter interval (l5j , both full reversible be- 
havior (Fig. and a complicated interrupted pattern 
(Fig. can be produced. In panels [2~^ and [2Jd the 
droplets were perpetually injected, such that the loop 
was not empty when changing the flow direction. Con- 
versely, in Fig. [2J; and [2ji, a finite train of N droplets 
was sent to and fro, with an empty loop in the mean- 
time. This latter case is hardly reversible and depends 
on the precise number of droplets. 

Let us underline that non-reversibility is deeply linked 
to the selection rule, and thus to the presence of junctions 
between transport channels. In a dual network with no 
such junctions, one expects the behavior to be reversible. 

Reversibility and Symmetry of the Design - This ob- 
servation leads us to a rather generic point regarding the 
kind of functions that a fore/aft symmetric device can 
achieve if the dynamics is reversible. A good support 
for this discussion is the recently proposed [7j] symmetric 
ladder network of Fig. [TJd. Two transport channels are 
connected by a few (narrower) bypass channels. The de- 
vice has been shown to perform the following function: 
Two droplets arriving with a given delay in the entrance 
channels leave the system with a much reduced distance, 
allowing synchronization between the two channels. Why 
is this remarkable? On one hand, (A) as the device is 
symmetric in design, if distances are reduced going one 
way, they should also be reduced going the reverse way. 
On the other hand, (B) the absence of junctions leading 
to "decision making steps" implies reversible dynamics, 
so that if the distance decreases one way, it should in- 
crease the other way: Such paradoxical statements are 
commonly used in low Reynolds number fluid dynamics 
to show that there can be no net effect. Here, experi- 
ments tell us otherwise. 

Running our model immediately resolves the paradox 
in a manner that is practically relevant: The function can 
be achieved only if the boundary conditions are not sym- 
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FIG. 3: Change of droplet misalignment in the ladder device. 
Only the boundary condition (a) leads to a contraction; fixed 
pressures on both ends (b) do not affect the offset, while it 
grows in (c). R d /R p =20, R hy /R p =5, L r /L p =5. 
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FIG. 4: Fraction of wrong (not strictly alternating) deci- 
sions in T-junctions with and without bypass (Fig. [it). The 
interval without errors around the average droplet distance 
(Axin) ~ 1-5L S corresponds to the experimental findings in 
Fig. 3 of Ref. The robustness of the behavior is chal- 

lenged by an asymmetry (L u /Le = 1.01) and a random per- 
turbation of the incoming distances (standard deviation is 5% 
of the average value). Rd/R s = 0.3. 



metric! For example, if we apply constant flow rates in 
the two entering arms and a fixed pressure at the two ex- 
iting nodes, we recover qualitatively the phenomenology 
described in Ref. p| with distance reduction (see case a 
in Fig. [3|. This reduction is not forbidden by the symme- 
try argument, as changing the sign of the flowrates and 
of the pressures yields a pattern of boundary conditions 
that is not the symmetric image of the initial situation, 
eliminating argument (A). In this new situation (case c 
in Fig. [3]) the ladder acts as a distance expander rather 
than a distance contractor, in accordance with (B). Con- 
versely, if we use fixed pressures at all end nodes, then 
the effect vanishes, as required by the symmetry argu- 
ment (A+B) above (case b in Fig. [3]). Our analysis thus 
clarifies the limits of use of this device as a synchronizing 
functional block in a network context. 

At this stage, we have demonstrated that our model 
reproduces qualitatively and quantitatively the perfor- 
mance of existing devices. It illuminates the essential 
role of the imposed driving conditions [l6j |. and it pro- 
vides a priori indication as to the dynamic behavior (pe- 
riod doubling, reversibility, etc.) as a function of design 
parameters (geometrical symmetry) and operational pa- 
rameters (droplets distance, etc.). 

Designing Devices: Optimization and Robustness - 
The speed of our numerical scheme permits fast opti- 
mization of identified designs. For example, we can scan 
the effect of different parameters of the ladder device 
(Fig. [3k.), such as the number of bypasses, their resis- 
tance, their distance, or the resistance of the droplets. 
The parameter Rd/R p has the strongest influence: the 
closer the bypasses, the better. Conversely, making the 
side arms of the transport channels (L r ) very long, has a 
negative effect on the contraction. Increasing the number 
of bypasses from 3 to 10, or making them 100 times less 
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resistant only weakly improves the contraction result. 

Beyond this optimization in an ideal world, an essential 
point for suggested designs to turn into efficient work- 
ing devices is their robustness to both fabrication errors 
and operational variations — virtues that are rarely quan- 
tified. We demonstrate the easy use of our tool to qualify 
the bypassed T-junction of Fig. [Ip, proposed in Ref. [l2|. 
The aim of this engineered junction is to provide reliable 
alternation of the channel chosen by incoming droplets. 
How it works is detailed in Ref. [12j - in short, the by- 
pass channel imposes similar pressures at both its ends. 
A droplet having chosen one route is then able to "block" 
its channel such that the next droplet favors the other one 
due to its larger flowrate. With our numerical scheme 
we immediately recover the same parameter window for 
which the experimental d esig n has been shown to be suc- 
cessful (cf. Fig. 3 in Ref. [l2(). 

We can systematically explore the robustness of the al- 
ternating behavior to imperfections of the device. These 
may consist of a "quenched" design asymmetry, such as a 
small systematic bias in the decision rule or a difference 
in arm lengths; or they may come as random fluctua- 
tions in the incoming droplet distance or in the pressures 
at the outlets. Figure H] shows the effect of a random 
Gaussian perturbation of the incoming droplet distances, 
combined with a slight asymmetry of the channel lengths 
behind the bypass. Both perturbations lead to the com- 
plete failure of an unbypassed junction, whereas even a 
thin bypass robustly yields a perfect alternation in a sig- 
nificant window of parameters. 

Playing with this analysis we find another remarkable 
property of the bypassed T-junction, namely noise re- 
duction in the droplet distances. The distances in the 
two outgoing sequences have a smaller standard devia- 
tion than the incoming train (data not shown). This 
effect operates for a single junction and can be amplified 
by using a cascade of several junctions, connected in such 
a way that the average hydrodynamic resistances are bal- 
anced at each stage of the cascade 17j. At each level of 
the cascade the noise is reduced because the relative stan- 
dard deviation of the distance between over-next droplets 
is only l/v2 times that between adjacent droplets. 

Conclusion - We have proposed a simple yet efficient 
fast numerical tool to analyze the traffic of droplets in 
"dual microfluidic" networks. Through an analysis of 
three recently proposed devices, we have shown the power 
of this tool for understanding fundamental issues, such as 
the role of the boundary conditions, of geometric sym- 
metry and of the occurrence of reversible behavior. We 
have left aside physical effects, in particular capillarity, 
that could add further functionalities, such as trapping 
in tapered channels, splitting and merging of droplets at 



junctions, etc. However, we believe, that the presented 
method will stimulate further experiments to explore the 
validity of its central ideas. After such a validation step 
it will provide a guide for the design of devices, and for 
exploration of new functionalities. In particular, it lends 
its ideas to robust passive solutions that are likely steps 
for droplet (or bubble) microfluidics to practically come 
up to the hopes it has raised in many fields. 

This work was supported by the French National Re- 
search Agency (ANR) via project Scan2, and we thank 
all members of this project for stimulating discussions. 
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